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(57) Abstract: A method of data estimation for a Time Division Duplex (TDD) Code Division Multiple Access (CDMA) system 
or any other system using an extended algorithm (EA) in preference to a truncated algorithm (TA). The EA avoids implementation 
errors by choice of proper extended matrices, and accepts the use of one piece of hardware. The EA also obviates loss of multiple 
signals in the tail part of each data field, and avoids errors due to transformation of a Toeplitz matrix to a circulant matrix. 
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[0001] EXTENDED ALGORITHM DATA ESTIMATOR 

[0002] FIELD OF INVENTION 

[0003] The invention generally relates to data estimation for wireless 
communication systems. More particularly, the invention relates to an extended 
algorithm (EA) for data estimation in such systems. 

[0004] BACKGROUND 

[0005] In some proposed wireless systems, data is wirelessly transmitted 
block by block with a separation interval between successive blocks. This 
property permits the application of joint detection (JD) in receivers to suppress 
the inter-symbol interference (ISI) and multiple-access interference (MAI). 
Single user detectors (SUDs) are used to estimate data of signals that go through 
a single downlink channel. An advantage of the SUD is that it can be 
implemented efficiently by fast Fourier transform (FFT), based on the rationale 
that a square Toeplitz matrix can be approximated as a circulant counterpart of 
the same size. 

[0006] When Toeplitz matrices are shortened along their longer dimensions 
to square matrices and replaced with their circulant counterparts, an 
approximation error is introduced. This error is most prevalent in the head and 
tail portions of the matrix. In many systems, the data associated with the head 
and tail portions carries system information required by the receivers, such as 
power control bits and the transport format combination indicator (TFCI) in the 
proposed third generation partnership project (3GPP) wideband code division 
multiple access (WCDMA) time division duplex (TDD) system. Itis desirable to 
enhance data estimation in such systems. 



[0007] SUMMARY 

[0008] The present invention provides a computationally efficient and 
accurate implementation of a data estimator for systems, such as Frequency 
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Division Duplex (FDD) or Time Division Duplex (TDD) Code Division Multiple 
Access (CDMA) systems. Described herein is a method of implementing a data 
estimator wherein no data of interest for estimation will be affected significantly 
by the circular approximation error. To achieve this, all of the square circulant 
matrices are extended. The advantages of the extended approach arise from two 
aspects: (1) the avoidance of loss of multipath signals in the tail part of each data 
field and (2) the avoidance of error due to the Toeplitz to circulant matrix 
transformation. As a result, longer discrete Fourier Transforms (DFTs) or FFTs 
are performed when implementing the EA. In order to minimize the required 
computations with DFT, the extended sizes should preferably be limi ted 
dynamically to their lower bounds according to specific data block length and 
channel delay spread. However, if prime factor algorithm (PFA) is used, 
increasing the FFT length usually does not result in the increase of 
computational complexity. It is to be noted that the computations can be 
minimized through choosing a proper FFT length in a certain range. In this case, 
fixed single-length FFTs ax^e desirable by considering the longest block length 
and delay spread. The single-length FFT with PFA makes it possible that 
different lengths of data block (burst types) are supported by one algorithm only. 
This simplifies the implementation further since one piece of hardware is needed 
to cope with the single algorithm. 

[0009] BRIEF DESCRIPTION OF THE DRAWING(S) 

[00 10] A more detailed understanding of the invention may be had from the 

following description of preferred embodiments given by way of example only and 

to be understood with reference to the accompanying drawing wherein: 

[0011] Figure 1 is a block diagram of a system used for implementing an 

EA with over-sampling in accordance with a preferred embodiment of the present 

invention; 

[00 12] Figures 2A, 2B and 2C, taken together, are a flow chart illustrating 
the method steps implemented by the EA of Figure 1. 
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[0013] Figure 3 is an illustration of raw block error rate (BER) of TFCI-1 
versus signal to noise (SNR) in Case 1 channel (TFCI represents transport format 
combination indicator); 

[0014] Figure 4 is an illustration of raw BER of TFCI-2 versus signal to 
noise ratio (SNR) per code in Case 1 channel; 

[0015] Figure 5 is an illustration of raw BER of all bits versus SNR per 
code in Case 1 channel; 

[0016] Figure 6 is an illustration of raw BER of TFCI-1 versus SNR per 
code in Case-2 channel; 

[0017] Figure 7 is an illustration of raw BER of TFCI-2 versus SNR per 
code in Case-2 channel; 

[00 18] Figure 8 is an illustration of all bits versus SNR per code in Case-2 
channel; and 

[0019] Figures 9A and 9B are receiver implementations using extended 
algorithm data detection. 

[0020] DETAILED DESCRIPTION OF THE PREFERRED 

EMBODIMENTS) 

[0021] The invention is generally applicable for data estimation in CDMA 
systems, such as 3GPP TDD mode and time division synchronous code division 
multiple access (TD-SCDMA). However, the following description, as an 
example, refers to a TDD CDMA system model and algorithm. In the described 
example, if a composite spread signal is sent from one transmitter site to one 
reception site, the received signal, r , is composite spread signal, 5, passed 
through a single channel, H_ . H is the channel response matrix. This process 
can be represented as r = ITs + n , where n is the noise vector. With W being the 
length of channel response, # takes the form per Equation 1. 
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. ft, . 

K-\ ■ • ■ 

A. 

Vi. 

Equation 1 

H is of size (L+W-l) by L. L denotes the number of chips in the time period of 
interest, such as a data field (block). The composite spread signal s can be 
expressed as s = C d , where the symbol vector d and the code matrix C take the 
form per Equation 2. 

d = (d l ,d 21 .., 1 d KNs )' 

Equation 2 

T denotes the transposition and C is per Equation 3. 

c=|c (1) ,c (2) ,...,c w J 

Equation 3 

Each C {k) is per Equation 4. 
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Equation 4 

[0022] Q, K, and N s (=L/Q) denote the spread factor (SF), the number of 

active codes, and the number of symbols carried on each channelization code, 
respectively. 

[0023] The SUD comprises two stages: (a) channel equalization and 0^) de- 
spreading. In the first stage, the composite spread signal s is estimated from 
r = Hs + n, preferably through a minimum mean squared error (MMSE) equalizer 
or a zero forcing solution. 
A MMSE equalizer is per Equation 5. 

S = [R H +<r 2 l] X H H r 
= [R H + <t 2 iYr h s + [r h + a 2 l] l H H n 

Equation 5 

A zero forcing solution is per Equation 6. 

I = Rh' 1 M" I > Equation 6 



C<*> = 



„<*> 
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I is the identity matrix, R H = H_ H # is a square Toeplitz matrix of size L, per 
Equation 7. 
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Equation 7 

' * * denoting the conjugate operation. In the second stage, a simple de-spreading 
process is performed to estimate the symbol sequence d, d, per Equation 8, 
d = C H s. 

Equation 8 

[0024] To implement Equation 5 efficiently, it is desirable to approximate 
the algorithm properly. To do this, first the Toeplitz matrix // is extended from 
the size of L+W-l by L to the size of Lm+W-1 by Lm and then the square matrix 
R H is extended from the size of L to Lm with Lm > L+W-l, while keeping the 
banded and Toeplitz structure of the matrices intact. The vector r is extended to 
the length Lm by zero padding, if the length of r is less than Lm. The vectors s 
and n are effectively automatically extended, due to the padding extension of the 
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other vector/matrix. The extended versions of the matrices and vector H,R H ,s, 
r and n are denoted as H_ E , R B> s E , r E and n E , respectively. R E is as follows, 

[0025] The last Lm-L elements of s E can all be deemed to be zero, which is 
essential in understanding the avoidance of the implementation error. With 
these notations, Equation 5 can be rewritten as Equation 9. 

l E = fes + <? 2 iYiIe Le> Equation 9 

r E is as follows r E = Ke£e + He ■ Because Equation 9 is the extended version of 

Equation 5, there should not be any differences between them if only the first L 
elements of S E are considered. The last W-l rows of the H_ E are cut to get a new 
square matrix of size Lm, denoted by H s . Suppose R ar and i£ Ctr represent the 
circular counterparts of R E and H_ s , respectively. From R E and # 5 , R ar and 
# Cir are constructed as per Equations 10 and 11. 







R WA 


0 


0 




R W-t 




K 


K 






R W-1 


0 












K 


'•• R , 






\ o 






Rw-i 








R t 




0 


0 






0 




r; 


Ro 






0 






0 


0 




k 






R «r-i 


0 






0 


'*• R w-\ 






"•• *i 




0 


0 






0 


R'wa 






R t 


R w-i 


0 






0 


0 


K-l 


r; 


R o •• 




R W-4 


R W-i 






0 


0 




R ; \ 


R, 












0' 






Ro 


Rt 






R W-l 






0 0 


R W-l 


r; 


Ro . 



Equation 10 
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. ft, 



Equation 11 



Rem and both of them are square matrices of size Lm. Because the last Lm- 
L elements of s E are zero, r E = H 5 s B + n E . Replacing R E , H H B and r E by R^, 
Kar and iL s i E +n E , respectively, in Equation 9, results in Equation 12. 



r , >, „ r t.1 „ Equation 12 

= War + VI] K"arIL s s E + [R ar +a 2 lfH H ar n E 

Assuming H s = H Cir -H A (H A is an error matrix), Equation 12 can be expressed 
as Equation 13. 



Equation 13 



Z = £La r H A s E is a column vector of length Lm and the error matrix H A is per 
Equation 14. 
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Equation 14 



The non-zero elements are located in a triangular area between the first W-l 
rows and the last W-l columns. Without the third term, Equation 13 is very 
similar to Equation 5 in function when the first L elements are considered. The 
vector y is evaluated. According to the structures of the matrices and # A , 

the matrix X_ = Kcu-Ka is a square matrix of size Lm with the structure per 
Equation 15. 



X = 



x ... x x 

X . X 



X X X 

X . X 

X X X 



Equation 15 



The non-zero elements, denoted by V, are located only in two areas: (1) a 
triangular one between the first W-l rows and the last W-l columns and (2) a 
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square one between the last W-l rows and columns. Because the last (Lm-L) 
elements of s E are zero, all elements of y = Ea r H A s E =Xs E are zero if 
Lm>L+W-L When Lm>L+W-l, Equation 16 results. 



Equation 16 is a good approximation of Equation 5, when taking only the first L 
estimates into account. The first portion of Equation 16 is referred to as the EA. 
Similarly, the head and tail parts of a data field are affected significantly by 
replacing Toeplitz matrices with their circulant counterparts in Equation 5 
directly without matrix extension. The implementation algorithm without 
matrix extension is called the truncated algorithm (TA). With TA, Equations 9 to 
15 are still valid, except that Lm=L. The reason consists of two aspects: First, 
when replacing H in Equation 5 with H Cir of size L, the length of the received 
signal vector r is preferably limited to be L. This results in the loss of multipath 
signals of the data in the tail part of a data field. Therefore, the estimates of the 
affected data become very poor. Second, when Lm=L, y is a column vector of 
length L, in which the first and last W-l elements are non-zero elements. 
Because B Cir = R Cir + a 2 1 is of banded structure along the diagonal, the inverse of 
B ar has approximately the same structure. Therefore, the relative large values 
of the column vector z = B'ly are located in the first and last W-l columns and 
affect the estimates of | both in the head and tail areas. For the second reason, 
the number of the affected estimates depends on the channel response length W. 
The larger the delay spread of a channel (W) is, the more the estimates are 
affected. 

[0026] Additionally, the implementation algorithm given by Equation 16 
can be extended to support over-sampling. With over-sampling, the effect of 
timing error is mitigated. Assuming sampling rate is M times of chip rate, M 



= fc* + a^Rol* + bar + ° 2 l\ X K 




Equation 16 
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received signal vectors, denoted by r £° for m=l,2, . . .,M, are available. However, 

the time interval between two successive samples in each r^ 0 is the chip 

duration. Similarly, there are also M sets of channel response, denoted by 

h im) =(/z 0w ,/4 m ,...,A w „ lOT )form=l,2,...,M. With these channel responses, atotalof 

2M circulant matrices H Grjtt and R ariM for m=l,2,...,M can be constructed. 

Accordingly, the implementation algorithm with over-sampling can be written as 
per Equation 17. 




Equation 17 



o 2 m is the noise variance corresponding to the mth input vector rf } . 

[0027] Before implementation, the value of Lm is determined. Because Lm 
is larger than L+W-l, Lm is chosen as per Equation 18. 

Lm = max {LRmax {W}+e, Equation 18 

max {.} is the maximum value of {.} and e is a positive integer used to let Lm a 
good length for FFT implementation. For example, in UTRA wideband TDD 
system (WTDD), max {L}=1104, max {W}=114. e is chosen to be equal to 14 to let 
Lm=1232. With this length, FFT can be performed very efficiently by PFA 
because 1232 can be factored as 1232=7x11x16. With complex input, the real 
multiplies and adds required by the 1232-point FFT are 8836 and 44228, 
respectively. From Equation 18, Lm depends on specific system design. 
However, the implementation approach described herein is applicable to any 
other TDD systems such as UTRA narrowband TDD system (TD-SCDMA). 
[0028] In the following, the preferred implementation procedure of 
Equation 17 is described as method steps under the assumption that the FFT 
length P equals the selected Lm. 

[0029] The first column g of the circulant matrix T f R r . + —a 2 1 1 is 

- al m m ) 

computed based on the estimated channel response and noise power, yielding 
Equation 19. 
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[0030] The circulant matrix YJ /U„ +^<rV\ ^ the FFT domain is 



WO 2004/023704 PCT/US2003/028163 

Equation 19 

decomposed , yielding Equation 20 

Equation 20 

and D" 1 are the P-point FFT and inverse FFT (IFFT) matrices defined as per 
Equation 21. 

P-l .j^ztn j p_j .2?ifa 

£V*= ' and D ?*=— » for k=0,l,.»P-l, 

n*0 ^ «=K> 



Equation 21 

A R is a diagonal matrix of size P, whose diagonal is A^is denoted as 

A* =diag(D p g). The relation between D' p and is D; 1 =(l/F)D*. 
[0031] The circulant matrix Kar^ is decomposed in the FFT domain, 

yielding Equation 22. 

Kor. m = d p h Hm D P , Equation 22 

A w ^ is a diagonal matrix of size P, whose diagonal is D p u m with 
l£m = Km > *W ' > Vm, .0. — ^ being the first column of H Qr>m . 
[0032] The received signal vector r (m) is reconstructed by zero padding to 

get the extended signal vector of length P. 

[0033] The composite spread signal vector 7 E is computed yielding 
Equation 23 or Equation 24 in the frequency domain. 

Equation 23 

-12- 
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° p1e = (|^^)^( D ^?0]/foi) and u = D?[DM, 

Equation 24 

The operators 0 and / denote the vector multiplication and division performed on 
element-by-element basis, respectively. The last P-L elements of J s are rounded 
off to get another vector | of length L. 

[0034] The composite spread signal s is despread, yielding? £ . 
[0035] Figure 1 is a block diagram of a system 100. For an over-sampled 

system, M sampled sequences are processed, r (1> . . . r iM) and h a) . . . ft <J ° . For a chip 
rate sampled sequence, only one sampled sequence is processed, r (1) and ft (1) . The 
system 100 receives the signals r (1> ... r (Af) at inputs 105i...l05 M (105) and 
receives the channel impulse response h m ...h m at inputs 110i ...110m (110). 
The received signals ... r (W) are zero padded in the tail by zero padding 
devices 115i. . .115m (115) until the length of each sequence achieves length Lm. 
The extended sequences after zero padding are denoted as r™ ... r£° which exits 
block 115 via outputs 120i...120m (120). The channel impulse responses 
h 0) ... h {M) are zero padded in the tail by zero padding devices 125i...125m (125) 
until the length of the extended sequence achieves length Lm. The extended 
sequences after zero padding are denoted as u } . . . u M which exits the zero padding 
devices 125 via outputs 130i...130m (130). DFT or FFT blocks 135i...135m (135) 
receive the outputs 120 from zero padding devices 115 and perform DFT or FFT 
on ...rj", F{r£ y )...F(r£ y ). DFT or FFT blocks 140i...140m (140) receive the 
outputs 130 from zero padding devices 125 and perform DFT or FFT on . . . u M , 
F{u x )...F(u M ). Conjugate devices 145i.. ,145m (145) conjugate F(u x ) ... F{u M ), 
F{u x )\.. F(u M ) * . Element-to-element multipliers 150i. ..150m (150) multiply the 
sequences FGff)...F(r«) and FQtJ ...F(u„Y , F{rfyFQ i O\..F(rf ) yF(u m )\ 
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[0036] All of the M sampled sequence results are added element-to- 

element by adder 175 with f^Fir^-FiuJ , M = 1,2,.. .,M. A channel 
correlation vector g is generated by a channel correlation generator 180 using 

M 

extended channel response sequences w t , with g-^g im} . 

[0037] Using a MMSE algorithm, a noise variance cr* is added to the first 

element of vector g {m) . Vector g <n,) is generated using u m . The i-th element of 
the vectors g (m) for the m-th sampled sequence is computed by: circulating the 

conjugate vector u m by downshifting i-1 elements and multiplying the shifted 
vector u m by the vector u m such that £ ( "\0=^M)^/b&«> w=l,2,...,Af . A DPT or 
PFT device 185 performs a DFT or FFT on channel correlation vector g , F(g) . 
Divider 190 divides element-by-element the output of adder 175 by the output of 

DFT or FFT device 185, such that . Inverse DFT or inverse 

FFT device 194 is performed on the output of divider 190, such that 

F~ l ) . The output of the inverse DFT or inverse FFT device 194 

Fig) 

is an estimate of the composite spread signal $ . Despreader 198 despreads the 
output of inverse DFT or inverse FFT device 194 to obtain the estimated data 
symbols d . 

[0038] Referring to Figures 2A, 2B and 2C, the procedure for performing an 
EA in accordance with a preferred embodiment of the present invention is 
described as follows: 

[0039] In step 205, system 100 receives signal r 0) at input 105 and receives 

the channel impulse response /i 0) at input 110. 
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[0040] In step 210, the received signal r (1) is zero padded in the tail by zero 

padding device 115 until the length of sequence achieves length Lm. The 
extended sequence after zero padding is denoted as which exits block 115 via 
output 120. 

[0041] In step 215, the channel impulse response h il) is zero padded in the 

tail by zero padding device 125 until the length of the extended sequence 
achieves length Lm. The extended sequence after zero padding is denoted as 
a, which exits the zero padding device 125 via output 130. 
[0042] In step 220, DFT or FFT block 135 receives the output 120 from 

zero padding device 115 and performs DFT or FFT on rf such that Firf) . DFT 
or FFT block 140 receives the output 130 from zero padding device 125 and 
performs DFT or FFT on w, , F(w j ) . 

[0043] In step 225, conjugate device 145 conjugates F(u x ) , F(u,)* . 

[0044] In step 230, element-to-element multiplier 150 multiplies the 

sequences F(rf) and F(m,)* producing F(r ( £ )-F(w 1 )*. 

[0045] In step 235, for over-sampling system with M sampled sequences, 

steps 210 to 230 are repeated for sampled sequences 2,. . .,M, F(r ( ^ } ) * F(u m )* , for m 
= 2,...,M. 

[0046] In step 240, all of the M sampled sequence results obtained in steps 

230 and 235 are added element-to-element by adder 175, £ F(r? } ) • F(u m )* , for m 
= 1,2,...,M. 

[0047] In step 245, a channel correlation vector £ is generated by a 

channel correlation generator 180 using extended channel response sequences 

M 

u, , such that g = ■ 

m=l 

[0048] In step 250, a DFT or FFT 185 performs DFT or FFT on channel 

correlation vector F(g). 
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[0049] In step 255, divider 190 divides element-by-element the result in 

£f<tP>-F<kJ- 



step 240 by the result in step 250, 



nt=l 



[0050] In step 260, an inverse DFT or inverse FFT 194 is performed on the 



result of step 255, F 



F(g) 



, producing the estimated composite 



spread signal, S . 

[0051] In step 265, despreader 198 despreads the output of step 260 to 

obtain the estimated data symbols d . 

[0052] In simulation, the described model is based on K = 12, and that 

codes are transmitted with equal code power and the effects of midamble signal 
on the second data field have been cancelled completely. Each code has a SF =16. 
A total 1104 chips in a data field are assumed (burst type 2 in WTDD). Because 
there are two data fields in a timeslot in WTDD, the last 8 bits (4 complex 
symbols) of the first data field, and the first 8 bits of the second data field are 
defined as TFCI-1 and TFCI-2, respectively. Two algorithms, truncated 
algorithm (TA) and extended algorithm (EA), are employed.. The raw bit error 
rates ofTFCI-1, TFCI-2 as well as all bits are evaluated for both EA and TA with 
chip rate sampling. 1000 timeslots are accumulated for each SNR point. The 
simulations are run over the WG4 Case-1 and 2 ch ann els. 
[0053] Figures 3 and 4 present the performance of TFCI-1 and TFCI-2 in 

WG4 Case-1 channel when EA and TA are used. A significant performance gap 
between EA and TA is found as shown in Figure 3. Because TFCI-1 locates in 
the tail part of the first data field, the performance degradation of the TFCI-1 
with TA results from two reasons: (1) loss of TFCI-1 multipath signals and (2) 
Toeplitz to circulant matrix replacement in the implementation because the 
channel response length W is small (W=4) for WG4 Case-1 channel. This 
conclusion is confirmed by the results shown in Figure 4. Because the 
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performance in Figure 4 is for the TFCI-2 which locates in the head part of the 
second data field, the most probable reason that affects its performance with TA 
must be the second one: matrix replacement. From Figure 4, it can be seen that 
the performance of TFCI-2 with EA and TA is nearly identical. This implies that 
the estimate error introduced to the TA through the matrix replacement is very 
limited because of the small value of W for WG4 Case-1 channel. For example, 
when W=4, only a quarter of the first symbol is affected because of the 
relationship SF=16. 

[0054] Figure 5 shows the raw BER of all bits when EA and TA are 

assumed. Comparing with EA, the major contributor to the loss of raw BER of all 
bits with TA is the TFCI-l in each timeslot. 

[0055] Figures 6 and 7 show the performance of TFCI-l and TFCI-2 in 
WG4 Case-2 channel when EA and TA are adopted. Case-2 channel differs from 
Case-1 channel by much larger delay spread (W=46) as well as stronger power of 
multipath signals. From Figure 6, the TFCM with TA is seen to be almost 
destroyed both due to the loss of its multipath signals and matrix replacement. 
In Figure 7, the performance of TFCI-2 with TA is still much poorer than that 
with EA, which is due to the matrix replacement only. In this case, W=46 and 
hence the first three symbols (six bits) will be affected significantly. Figure 8 
presents the raw BER of all bits for EA and TA in Case-2 channel. 
[0056] When a truncated algorithm (TA) is used in the implementation, the 
data in the head and tail parts of a data field are affected significantly by two 
aspects: the lost information of multipath signals due to cutting the channel 
response matrix to be square and the error due to replacing Toeplitz matrices 
with circulant ones. To overcome this problem, an EA is used. The EA avoids 
implementation errors by choosing the size of the extended matrices properly. To 
implement the EA with DFT, a dynamic-length EA is desirable, while for an EA 
with FFT of PFA, a fixed-length approach is more appropriate. In the fixed- 
length approach, the computation complexity can be minimi zed through choosing 
a proper PFA length in a certain range. The fixed-length EA makes different 
data block lengths (burst types) be supported by one algorithm only. The fixed 
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length EA simplifies the implementation further since one piece of hardware is 
needed to cope with the single algorithm. Simulation results show that the 
performance of the EA is much better than that of the TA, especially for the data 
in the head and tail parts of data fields. 

[0057] The invention can be implemented at a base station or wireless 
transmit/receive unit (WTRU). Hereafter, a wireless transmit/receive unit 
(WTRU) includes but is not limited to a user equipment, mobile station, fixed or 
mobile subscriber unit, pager, or any other type of device capable of operating in 
a wireless environment. When referred to hereafter, a base station includes but 
is not limited to a base station, Node-B, site controller, access point or other 
interfacing device in a wireless environment. 

[0058] Figures 9A and 9B are receiver implementations using extended 
algorithm data detection. Referring to Figure 9A, radio frequency (RF) signals 
are received by an antenna 300. A sampling device 305 produces a chip rate 
received vector r. A channel estimation device 325 determines a channel impulse 
response h for the received vector. A single user detection device 310 uses the 
received vector r and the channel impulse response h to estimate the data vector 
d using the extended algorithm. The received vector r is processed by a channel 
equalizer 315 using the channel impulse response h to determine a spread vector 
s. A despreader 320 using the transmission codes C despreads the spread vector 
s to estimate the data vector d. 

[0059] Referring to Figure 9B, RF signals are received by an antenna 300. 
A sampling device 305 samples the received signal at a multiple M of the chip 
rate, producing M received vector sequences a. . -Cm. A channel estimation device 
325 determines a channel impulse response Ai...Am corresponding to each 
received vector a • - .£m. A single user detection device 3 10 uses the received vector 
sequences a. . .r m and the channel impulse responses & . . .h m to estimate the data 
vector d using the extended algorithm. The received vectors o- • -Cm are processed 
by a channel equalizer 315 using the channel impulse responses hh..h m to 
determine a spread vector S. A despreader 320 using the transmission codes C 
despreads the spread vector S to estimate the data vector d. 
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[0060] While the present invention has been described in terms of the 
preferred embodiment, other variations which are within the scope of the 
invention as outlined in the claims below will be apparent to those sMlled in the 
art. 

* * * 
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CLAIMS 

What is claimed is: 

1. A method for recovering data from a plurality of signals received 
in a shared spectrum, the plurality of signals experiencing a similar channel 
response, the method comprising: 

sampling a composite signal including the plurality of received signals, 
producing a received vector; 

estimating a channel response of the composite signal; 
extending the received vector; 
extending the channel response; 

channel equalizing the received vector using the extended channel 
response, producing a spread vector; and 

despreading the spread vector to produce data of the plurality of 
signals. 

2. The method of claim 1 wherein a time interval between two 
successive samples in each extended received vector is the chip duration. 

3. The method of claim 1 wherein a time interval between two 
successive samples in each extended received vector is a fraction of the chip 
duration. 

4. The method of claim 1 farther comprising: 

computing a first column of a circulant matrix based on estimated 
channel response and noise power; 

decomposing a received vector circulant matrix in a fast Fourier 
transform (FFT) domain; 

decomposing a channel response circulant matrix in the fast Fourier 
transform (FFT) domain; 
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reconstructing the received signal vector resulting in an extended signal 
vector; 

computing the composite spread signal vector; and 
despreading the composite spread signal. 

5. A base station including a communications receiver, the receiver 
comprising: 

an antenna for receiving radio frequency (RF) signals; 
a sampling device coupled to the antenna for producing a chip rate 
received vector; 

a channel estimation device coupled to the sampling device for 
determining a channel impulse response for the received vector; and 

a single user detector (SUD) coupled to the sampling device and the 
channel estimation device for estimating a data vector using an extended 
algorithm which extends the received vector and the channel impulse 
response. 

6. The base station of claim 5, wherein the SUD comprises: 
a channel equalizer for using the channel impulse response to 

determine a spread vector; and 

a despreader coupled to the channel equalizer for despreading the 
spread vector to estimate the data vector. 

7. A wireless transmit/receive unit (WTRU) including a 
communications receiver, the receiver comprising: 

an antenna for receiving radio frequency (RF) signals; 
a sampling device coupled to the antenna for producing a chip rate 
received vector; 

a channel estimation device coupled to the sampling device for 
determining a channel impulse response for the received vector; and 
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a single user detector (SUD) coupled to the sampling device and the 
channel estimation device for estimating a data vector using an extended 
algorithm which extends the received vector and the channel impulse 
response. 

8. The WTRU of claim 7, wherein the SUD comprises: 

a channel equalizer for using the channel impulse response to 
determine a spread vector; and 

a despreader coupled to the channel equalizer for despreading the 
spread vector to estimate the data vector. 

9. A base station including a communications receiver, the receiver 
comprising: 

an antenna for receiving radio frequency (RF) signals; 

a sampling device coupled to the antenna for sampling the received 
signals at a multiple M of the chip rate, producing M received vector 
sequences; 

a channel estimation device coupled to the sampling device for 
determining a channel impulse response for each received vector; and 

a single user detector (SUD) coupled to the sampling device and the 
channel estimation device for estimating a data vector using an extended 
algorithm which extends the received vector and the channel impulse 
response. 

10. The base station of claim 9, wherein the SUD comprises: 

a channel equalizer for using a channel impulse response to determine 
a spread vector; and 

a despreader coupled to the channel equalizer for despreading the 
spread vector to estimate the data vector using transmission codes in the 
received signals. 
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11. A wireless transmit/receive unit (WTRU) including a 
communications receiver, the receiver comprising: 

an antenna for receiving radio frequency (RF) signals; 

a sampling device coupled to the antenna for sampling the received 
signals at a multiple M of the chip rate, producing M received vector 
sequences; 

a channel estimation device coupled to the sampling device for 
determining a channel impulse response for each received vector; and 

a single user detector (SUD) coupled to the sampling device and the 
channel estimation device for estimating a data vector using an extended 
algorithm which extends the received vector and the channel impulse 
response. 

12. The WTRU of claim 11, wherein the SUD comprises: 

a channel equalizer for using a channel impulse response to determine 
a spread vector; and 

a despreader coupled to the channel equalizer for despreading the 
spread vector to estimate the data vector using transmission codes in the 
received signals. 

13. A single user detector (SUD), comprising: 

(a) a channel equalization stage, wherein a composite spread signal is 
estimated using a minimum mean squared error (MMSE) equalizer; and 

(b) a despreading stage for estimating symbol sequences detected by the 

SUD. 

14. A communications system comprising: 
a base station; and 

a wireless transmit/receive unit (WTRU) in communication with the 
base station, wherein the WTRU comprises: 

an antenna for receiving radio frequency (RF) signals; 
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a sampling device coupled to the antenna for producing a chip rate 
received vector; 

a channel estimation device coupled to the sampling device for 
determining a channel impulse response for the received vector; and 

a single user detector (SUD) coupled to the sampling device and the 
channel estimation device for estimating a data vector using an extended 
algorithm. 

15. The communications system of claim 14, wherein the SUD 
comprises: 

a channel equalizer for using the channel impulse response to 
determine a spread vector; and 

a despreader coupled to the channel equalizer for despreading the 
spread vector to estimate the data vector using transmission codes in the 
received signals. 

16. A communications system comprising: 
a wireless transmit/receive unit (WTRU); and 

a base station in communication with the WTRU, wherein the base 
station comprises: 

an antenna for receiving radio frequency (RF) signals; 

a sampling device coupled to the antenna for producing a chip rate 
received vector; 

a channel estimation device coupled to the sampling device for 
dete rmini ng a channel impulse response for the received vector; and 

a single user detector (SUD) coupled to the sampling device and the 
channel estimation device for estimating a data vector using an extended 
algorithm. 

17. The communications system of claim 16, wherein the SUD 
comprises: 
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a channel equalizer for using the channel impulse response to 
determine a spread vector; and 

a despreader coupled to the channel equalizer for despreading the 
spread vector to estimate the data vector using transmission codes in the 
received signals. 

18. In a wireless communication system, a method for performing an 
extended algorithm (EA) with over-sampling, the method comprising: 

t (a) the system receiving a signal r (1) at a first input and a channel impulse 
response h 0) at a second input; 

(b) zero padding the received signal r a) in the tail until the length of 
sequence achieves length Lm and denoting the extended sequence after zero 
padding as rf ; 

(c) zero padding the channel impulse response in the tail until the 
length of the extended sequence achieves length Lm and denoting the extended 
sequence after zero padding as u x ; 

(d) performing a discrete Fourier Transform (DFT) or fast Fourier 
transform (FFT) on rf such that F(rf); 

(e) performing DFT or FFT on k, such that F(u, ) ; 

(f) conjugating F(u x ) such that F(u x y ; 

(g) multiplying the sequences F(r ( j?) and F(u x )* such that F(r£ ) VF(tt l )\ 
wherein for M sampled sequences, steps (b) - (g) are repeated for sampled 
sequences 2.....M such that F(r% } )- F(u m ? , m = 2,..., M. 

19. The method of claim 18, wherein all of the M sampled sequence 
results obtained in steps (b) - (g) are added element-to-element such that 

^FCr^J-FfeJ* ,M=1,2,...,M. 
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20. The method of claim 19 further comprising: 

(h) generating a channel correlation vector g using extended channel 

M 

response sequences w, such that g = g (w> ; 

m=l 

(i) performing DFT or FFT on channel correlation vector g such that 
(j) dividing element-by-element the result in step (g) by the result in step 

M 

(i) such that ■ 2s! ; 

F(g) 

(k) performing an inverse DFT or inverse FFT on the result of step (j) 

such that F~\m );and 

(1) despreading the result of step (k) to obtain the estimated data 

A 

symbols d . 



21. A wireless communication system for performing an extended 
algorithm (EA) with over-sampling, the system comprising: 

(a) means for receiving a signal r (1) at a first input and a channel impulse 
response h il) at a second input; 

(b) means for zero padding the received signal r (1> in the tail until the 
length of sequence achieves length Lm and denoting the extended sequence after 
zero padding as rf ; 

(c) means for zero padding the channel impulse response h 0) in the tail 
until the length of the extended sequence achieves length Lm and denoting the 
extended sequence after zero padding as w, ; 
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(d) means for performing a discrete Fourier Transform (DFT) or fast 
Fourier transform (FIT) on rf such that F(rf); 

(e) means for performing DFT or FFT on u x such that F{u x ) ; 

(f) means for conjugating F{u x ) such that F(u x )* ; 

(g) means for multiplying the sequences F(rf) and FQtJ such that 
Fte.™)* F(u x )* , wherein for M sampled sequences, steps (b) - (g) are repeated for 
sampled sequences 2,...,M such that F(r { ^)'F(u m )\ m = 2,..., M. 

22. The system of claim 21, wherein all of the M sampled sequence 
results are added element-to-element such that F(r£°) ■ F(u m )* , M = 1,2, . . .,M. 

23. The system of claim 22 further comprising: 

(h) means for generating a channel correlation vector g using extended 

M 

channel response sequences w, ,...,u M such that g =^g im) ; 

(i) means for performing DFT or FFT on channel correlation vector g such 
that F(g); 

(j) means for dividing element-by-element the result in step (g) by the 

result in step (i) such that ^ : 

F{g) 

(k) means for performing an inverse DFT or inverse FFT on the result of 

f,F(r™)-F(Hj 
step (j) such that F'\^ ) ; and 

(1) means for despreading the result of step (k) to obtain the estimated data 

A 

symbols d. 



-27- 



WO 2004/023704 



PCT/US2003/028163 



(5 




WO 2004/023704 



PCT/US2003/028163 



FIG, 2 

FIG.2A 
FIG. 2B 
FIG. 2C 



2/8 




RECEIVE SIGNAL rO) AND CHANNEL 
IMPULSE RESPONSE 



FIG.2A 



ZERO PAD THE RECEIVED SIGNAL fO) 
IN THE TAIL UNTIL THE LENGTH OF SEQUENCE 
ACHIEVES LM. DENOTE THE EXTENDED SEQUENCE 
AFTER ZERO PADDING AS rf ^ 



ZERO PAD THE CHANNEL IMPULSE RESPONSE h( 1 ) 

IN THE TAIL UNTIL THE LENGTH OF SEQUENCE 
ACHIEVES LM. DENOTE THE EXTENDED SEQUENCE 
AFTER ZERO PADDING AS IJ1. 



> 

PERFORM DFT OR FF 

WERE OBTAINED IN 
RESPECTIVELY SUCH 


r 

r ON r0) AND U-| THAT 
STEPS 210 AND 215 
THAT F(r0)) AND F(U-j). 


> 


f 



PERFORM CONJUGATE ON F ( yi ) SUCH THA T F(U 1 )* 



•220 



225 



WO 2004/023704 



PCT/US2003/028163 



3/8 



PERFORM ELEMENT-TO-ELEMENT 
MULTIPLICATION OF SEQUENCE 
F(rd)) AND ?(Utf SUCH THAT F([W) ■ F(m)* 



FOR OVER-SAMPLING SYSTEM WITH M SAMPLED 
SEQUENCES, REPEAT STEPS 21 0 TO 230 FOR 

SAMPLED SEQUENCES 2, M SUCH THAT 

F(K m ))-F(u m )*,m=2,... I M. 



PERFORM ELEMENT-TO-ELEMENT ADDITION 

OF RESULTS OBTAINED IN STEPS 230 AND 235 

FOR ALL M SAMPLED SEQUENCES SUCH THAT 
M 

2F(r_N). F (U m )\M=1,2 M. 

m=l E 



GENERATE CHANNEL CORRELATION VECTOR g 

USING EXTENDED CHANNEL RESPONSE 

M 

SEQUENCES ^ u M SUCH THAT g = £ g( m ) . 

" m=i - 



230 



■235 



240 



■245 



FIG. 2B 



WO 2004/023704 



PCT/US2003/028163 



4/8 



> 


r 


PERFORM DFTOR FFTON CHANNEL CORRELATION 
VECTOR g SUCH THAT F(g). 




> 


f 


DIVIDE ELEMENT-BY-ELEMENT THE RESULT IN 

STEP 240 BY THE RESULT IN STEP 250 SUCH THAT 
M 

2F(r_M) • F (u m r 
m=i E 

F(g) 




> 


r 


PERFORM INVERSE DFTOR INVERSE FFTON 

THE RESULT IN STEP 255 SUCH THAT 
M 

2F(rN). f (u m r 
pl( m=l E 

{ F(g) ' 




> 


r 


DESPREAD THE RESULT IN STEP 260 TO OBTAIN 
THE ESTIMATED DATA SYMBOLS 3. 





250 



255 



260 



265 



END 



FIG. 2C 



WO 2004/023704 



PCT/US2003/028163 



5/8 



rawBER 
of TFCl-1 




SNR PER CODE [dB] 

FIG. 3 

. WG4 CASE 1, 12 CHANNELIZATION CODES, SF = 16 EACH 



rawBER 
0fTFCI-2 10' 



\l 1 "~ 1 1 1 1 I 






8^ i i i i i 


I 


^^K^ 1111 
_ l ^s.- 1 1 -1- J 


1 1 
1 

) 


i ^"tS. i » 1 

1 1 I t 

1 1 1 1 » 

] 1 1 1 


1 


£ 1 H + 




I 1 1 1 X 


t 


lllll x. 


~B- TRUNCATED ALGORITHM 
EXTENDED ALGORITHM^ 


III! 
ill) 
llli 



-4 -2 0 2 4 6 8 10 12 



SNR PER CODE [dB] 

FIG. 4 



WO 2004/023704 



PCT/US2003/028163 




10° 



WG4 CASE 2, 12 CHANNELIZATION CODES, SF = 16 EACH 



-1 



10 



rawBER -jq-2 
Of TFC1-1 




10- 3 ; = , = = =: = = = = J=: = -_ = = 



io- 



-B- TRUNCATED ALGORITHM 
-fc- EXTENDED ALGORITHM 



2 4 6 
SNR PER CODE [dB] 



10 



12 



FIG. 6 



WO 2004/023704 



PCT/US2003/028163 



7/8 




SNR PER CODE [dB] 

FIG. 7 



rawBER 
of all bits 



10 



-1 



WG4 CASE 2, 12 CHANNELIZATION CODES, SF = 16 EACH 



10"' 



10" 















r 3 n _ _ ~ _ _ 




1 




i t 


I 

1 
1 


nEEr 




i 

. i 


i _ j 




i i 
i \ 








1 


l t 


i t 


1 

1 


i I 
1 i 


i X. i 










1 

1 






""i 

--1 


\. "i - -* -Sr. 

i _ _ J _ _ _ 3 


. -B- TRUNCATED ALGORITHM 
EXTENDED ALGORITHM 


i 




i i 
i i 

i »- ..... 



10^ = = = = : 



2 4 6 
SNR PER CODE [dB] 

FIG. 8 



10 



12 



WO 2004/023704 PCT/US2003/028163 



8/8 



300 



L 



305 



SINGLE USER DETECTION DEVICE 310 

C 



315 



SAMPLING 
DEVICE 


r 




CHANNEL 
EQUALIZER 


S 






> 



320 



DESPREADER 



CHANNEL 
ESTIMATION 
DEVICE 



325 



-> d 



FIG. 9 A 




300 



r 



305 



SAMPLING 
DEVICE 



11 



IM 



SINGLE USER DETECTION DEVICE 310 

C 



315 



CHANNEL 
EQUALIZER 



• • • # 



hi 



LL 



320 



DESPREADER 



>| CHANNEL 
ESTIMATION 
DEVICE 



•325 



> d 



F/G. 98 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE^) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: . . 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



